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Abstract. We present a theoretical study of the bandstructure and Landau levels 

in bilayer graphene at low energies in the presence of a transverse magnetic field and 

Rashba spin-orbit interaction in the regime of negligible trigonal distortion. Within 

an effective low energy approach (Lowdin partitioning theory) we derive an effective 

Hamiltonian for bilayer graphene that incorporates the influence of the Zeeman effect, 

1 . , the Rashba spin-orbit interaction, and inclusively, the role of the intrinsic spin-orbit 

^ | interaction on the same footing. Particular attention is spent to the energy spectrum 

Q ■ and Landau levels. Our modeling unveil the strong influence of the Rashba coupling 

Xr in the spin-splitting of the electron and hole bands. Graphene bilayers with 
weak Rashba spin-orbit interaction show a spin-splitting linear in momentum and 
proportional to Ajj, but scales inversely proportional to the interlayer hopping energy 
7i. However, at robust spin-orbit coupling A^, the energy spectrum shows a strong 
Q*\ ' warping behavior near the Dirac points. We find the bias-induced gap in bilayer 

graphene to be decreasing with increasing Rashba coupling, a behavior resembling a 
topological insulator transition. We further predict an unexpected assymetric spin- 
splitting and crossings of the Landau levels due to the interplay between the Rashba 
interaction and the external bias voltage. Our results are of relevance for interpreting 



magnetotransport and infrared cyclotron resonance measurements, including also 
situations of comparatively weak spin-orbit coupling. 



X 

H 

d ■ PACS numbers: 73.22.Pr, 75.70.Tj, 71.70.Di, 71.70.Ej 
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1. Introduction 

Graphene and its bilayer (BLG) posses very distinctive physical properties. [U E] 
Neglecting interactions, the low energy quasiparticles in pristine single and double layer 
graphene obey, respectively, linear and quadratic dispersion laws at the K(K') symmetry 
points. jH El E] In the presence of a quantizing magnetic field B > 0, the relativistic 
massless Dirac fermions of graphene exhibits Landau levels (LLs) non-equidistant in 
energy, E n ex \J\n\B, with n — 0, ±1, ±2, ..., the LL index. [6]The latter gives rise to the 
half integer quantum Hall effect at room temperature, [7] and to the fractional quantum 
Hall effect at high magnetic fields in the presence of many body effects. [HI E] On the 
other hand, LLs for BLG show a rather intricate index sequence instead; with a roughly 
linear B- field dependence for low LLs, and a \[B for high LLs. pm [HI 02] At low energies, 
the LLs in unbiased BLG follows the sequence E n oc a/M(M + 1) for n > 1 with a 
double degenerate zero-energy level, E = for n = 0. [131 US] Experimentally, 
the LLs dipole-allowed transition energies in single layer and BLG have been studied 
in detail by cyclotron resonance. [161 El HE] Most recently, phonon-mediated inter LL 
excitations have been explored by magneto- Raman scattering experiments. [19] 

The particularly high interest in graphene spin physics is strongly motivated by 
its expected prospects in nanoelectronics and spin-based devices for spintronics. [20] In 
this realm, the role of the spin-orbit interaction (SOI) effects in graphene sheets is a 
phenomenon under intense scrutiny. Two types of SOI in graphene have been identified; 
(«) the induced by carbon intra-atomic SOI (intrinsic-SOl) , and (ii), the SOI coming 
from the breaking of the space inversion symmetry of the hexagonal lattice (extrinsic- 
SOI), this is the so-called Rashba-SOI. The latter can have different origins, among 
those is the presence of a substrate, buckling, ad-atoms, or external electric fields. 

The magnitude of the excitation gap rj of the intrinsic-SOI in single layer is 
predicted to be very small (0.86-50 /xeV). [2D 1221 1231 ES ESI EEJ [27] Likewise, estimates of 
the Rashba coupling X R , leads to small values (few tens of /xeV) at typical electric fields 
(~ 0.16 mV/nm). [23 l l24 l 125] However recent spin-resolved photoemission experiments 
in graphene/Au/Ni(lll) showed enhancements of the Rashba coupling as large as 13 
meV. [28] Induced distortions by neutral impurities (ad-atoms) [291 SO] and the interplay 
of buckling with Rashba-SOI also yield significant enhancement of the spin-splittings 
(up to 40meV).|31j The role of the impurities and lattice deformations seems to be 
crucial for the observed long spin relaxation times (up to 2 ns in BLG[32]) linked to SOI 
effects[33l [31], as well as for a non-monotonic dependence on bias voltage of the spin 
relaxation times in BLG due to the D'yakonov-Perel' spin-precession mechanism. [3.5] 
Large spin-splittings (~ 0.22 eV) of the graphene 7r— states attributed to Rashba-SOI 
has been reported also in epitaxial graphene on a Ni substrate. [36] In theory, it has 
been shown that the Rashba-SOI can induce significant changes in the bandstructure of 
graphene [3?1 [38] as well as in BLG. [39] 

In graphene monolayers with Rashba coupling, the Landau levels are described by 
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Eil± = M\/^a> ( in units of X R ), with[ 



£<& = (2n - l)f 2 + -\l±y/l + A(2n - l)r 2 + 4f M (1) 

for n > 2, where f 1 = ^up/Z^IA^, t^ is the Fermi velocity (~ 10 6 ms~ 1 ), Ib = y/ch/eB 
is the magnetic length, h the Plank's constant over 2n, c is the light velocity in vacuum, 
and — e is the electron charge. Here /i = ± stands for the electron/hole LL branch. The 
lowest n = is given by Eq = whiles the n = 1 level gives rise to three modes. A zero 
mode E{ 0) = 0, in addition to its two nondegenerate modes at Ex^ = [iy 1 + 2r 2 . Exact 
solutions for the LLs in monolayer graphene under the influence of a Zeeman field and 
spin-orbit interactions has been also reported recently by De Martino et al. . |40] 

In this paper we show within low energy effective theory that for biased BLG in 
which the Rashba effect is the dominant SOI, its LLs must follow E { ^ ± = /Ua/£1 2) ±, with 

S^l = U 2 + | (r 2 + 2n u? ± V4u; 4 + 4nw 2 r 2 + r 4 ) (2) 

for n > 2, with U the interlayer bias energy, Y = 2y/2\ R VFh/'jilB, and to = 2v 2 7 h 2 /'jil B , 
being 71 inter-layer hopping energy. As it occurs in graphene, in BLG with Rashba-SOI 

we obtain three modes for n = 1, namely, the nondegenerate E{J + = fi\/U 2 + T 2 + 2u 2 

(2) 
and E\_ = —U, whereas for n = its eigenvalue also vanishes in the absence of gating, 

Eq} = U. Eq.(T5]) comprises one of the main results of this contribution. 

The aim of this study is to investigate the energy spectrum and the Landau levels 
in bilayer graphene under the influence of sizable spin-orbit interactions (SOI) of both, 
intrinsic and Rashba types. An effective low energy Hamiltonian for bilayer graphene 
that includes both types of SOI and Zeeman effect is derived within Lowdin partitioning 
theory. Whiles the Rashba SOI in single layer graphene is known to modify its otherwise 
conic spectrum, to a spectrum that includes two zero gap bands and two gapped 
branches of width 2X R (with parabolic shape, similar to unbiased BLG);[37J here in 
contrast, the (unbiased) BLG with Rashba-SOI shows a spin-splitting which is linear in 
momentum and proportional to \r, but inversely proportional to the interlayer hopping 
parameter 71. We predict also a strong influence of the Rashba spin-orbit interaction 
in the warping of the low energy bandstructure of biased bilayer at comparatively 
weak spin-orbit coupling \r. It is found that the bias-induced gap in bilayer graphene 
decreases as the Rashba strength coupling is increased. It is further predicted the rise 
of an unexpected asymmetric spin-splitting of the Landau levels due to the interplay 
among the Rashba coupling and the external bias voltage. 

The remainder of the paper is organized as follows. In Sec. 2 we outline the 
model and derivation for the low energy BLG effective Hamiltonian. The Landau level 
spectrum in the presence of Rashba-SOI is discussed in Sec. 3. The band spectrum 
properties and the Landau levels of BLG with Rashba-SOI are analyzed in detail in 
Sec. 4 and Sec. 5. A summary of our results is given in Sec. 6 . Additionally, we provide 
three appendixes. In Appendix A we outline the derivation (Lowdin partitioning theory) 
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of the low energy Hamiltonian in the presence of intrinsic and Rashba types of SOI, as 
well as the Zeeman effect. In Appendix B a detailed description of the eigenvalues of 
the low energy Hamiltonian is given, and finally in Appendix C the Landau levels for 
BLG with Rashba-SOI are derived. 



2. Low energy effective Hamiltonian for bilayer graphene 

Here we focus in the derivation of the low energy effective Hamiltonian for BLG 
with SOI in the presence of magnetic field that eventually leads to Eq.(j2]). We 
start by considering a pile of two graphene layers (BLG) in which the sites A 2 of 
the upper layer 2 lies directly on top of the B\ sites of the bottom layer 1 (AB- 
Bernal stacking). At the vicinity of the Dirac K symmetry point, the effective non- 
interacting bilayer graphene Hamiltonian H , written in terms of the spin- dependent 
basis \^ K ) = {i)A 1 +,i>A u ,i>B lf ,i>B u ,i>A n ,i>A n ,i>B2 r ,'il>B n }> nas the 8 x 8 matrix form 



( ±u 




H n 



where n = n x + iir v , with 7r 



H 4 



\ 





77T 





'j- 





±u 



77T 
(7T,,, 7T, 



77T 



±u 





77T 





\ 

t 



(3) 



xi "yj 



±u J 

is the canonical momentum, 



p + eA/c 

and A is the vector potential. Here 7 = Vp = , y a\/3/2h, with 7„ ~ 3.1 eV (intra- 
layer hopping energy) and a = 0.246 nm is the lattice parameter. pQ The electrostatic 
potential ±U of the bottom/upper layer is gate voltage adjustable and opens a gap of 
2U in the spectrum. [2] The dominant interlayer interaction in BLG is described to first 
approximation by the term 

/ — v^ 71 

-v^ 

Q IIaTT* 

V 3 7T 



W 



\ 



V 



v 3 ir 




-W47T' 





7i 




(4) 



-v^tc' J 



where 71 is the nearest neighbor (interlayer) hopping energy (~ 0.l7 o ). The terms 
proportional to the velocities v 3 = j 3 a\/ r 3/2h and t> 4 = j 4: a\/ r 3/2h arise due to 
second nearest neighbor (interlayer) hopping processes associated with 73 and 74 tight- 
binding parameters, respectively. |15[ |4"T] The former produces a trigonal warping whose 
characteristic energy is E 3 = ^m*v 3 2 = ■y 1 ( r y 3 /2 r y ) 2 ~ 1 meV, being m* = 7 1 /2uf, the 
electron effective mass [15], whiles the latter has yet a smaller characteristic energy, 
E A = \m*v\ = 7i(74/27 ) 2 ~ 0.2 meV. Since typically E 3A < 71 in BLG, the 
contributions in Eq.fllJ) coming from the terms with velocities v 3 and V4 are negligible. 
However we should notice that, when considering Rashba spin-orbit interaction, there 
might be situations in which v 3 and V4 may be important, particularly at very weak 
Rashba-strengths of the order of E 3 ^. We would like to mention nevertheless, that 
there is experimental (and theoretical) evidence of Rashba-SOI spin-splittings of one 
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order of magnitude [2H1 1221 HD] and even larger [35] than the typical distortion energies 
£3,4. In any case, its inclusion in the model can be readily incorporated in the general 
derivation of the reduced effective Hamiltonian via the Lowdin perturbation theory 
described in Appendix A. The regime in which there is a possible interplay among the 
1*3,4 terms and Rashba SOI is out of the scope of the present analytical study and for 
sake of clarity and simplicity these trigonal terms will be disregarded hereafter. 

When taking into account explicitly the presence of intrinsic- SOI, Rashba-SOI and 
Zeeman splitting, the total eight-band effective Hamiltonian will read 



Hr = Ho + Hr + Hi + Hz- 



(5) 



The second term to the right in Eq.(jSJ) arises due the influence of an effective electric 
field perpendicular to the BLG plane producing a Rashba type of SOI. The the leading 
contribution to the Rashba-SOI Hamiltonian, Hr, is modeled as follows: 
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i\ R a 
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-iX R a+ 






\ 



/ 



(6) 



here \r parameterize the strength of the intra-layer Rashba-SOI, as in monolayer 
graphene, with <J± = \ ((J x ±i& y ), being (a x ,cr y ) the usual 2x2 Pauli spin matrices. 
The intensity of the Rashba-SOI can be sizable (A# ~ 10 meV) due for instance to 
the presence of a metallic substrate. [28] Within tight-binding theory, it is understood 
that the Rashba-SOI arises because of the effective nearest-neighbor hopping of two p z 
orbitals with opposite spins under the presence of an applied transverse electric field. [26] 

Recently it has been predicted that Xr can be even larger (of a few tens of meV) 
due to buckling effects in conjunction with external electric fields. [31] Furthermore by 
varying the electric field the Rashba parameter can be tuned. A possible inter-layer 
Rashba spin-orbit coupling of strength A^ can in principle be present in BLG as well, 
however such contributions will be ignored here because of its predicted slight influence 
on the energy bands for Xr/'Jo ^ 0.3. |39j 

Additionally, in the same basis set above, the intrinsic-SOI Hamiltonian for BLG 
(third term to the r.h.s. in Eq.(jS]) ) shall follows the 8x8 matrix form 

\ 



T]S Z 

—7] S 



Hi 



I TjS z 
















(7) 



/ 



with i] the intrinsic SOI constant and s z is the spin operator along the z-axis, 
perpendicular to the BLG plane. The intrinsic SOI is a second order tunneling process 
(within tight-binding theory), which involves next-nearest-neighbor hopping events of p z 
electrons of a given spin. As mentioned in the introduction, the value of the excitation 
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energy r] has been inferred to be very small in monolayer graphene in both K(K') 
symmetry points (< 50/ieV), even if one consider interactions up to first order by 
including the unoccupied d and higher orbit als. j25l |26] Interestingly, in BLG, taking 
into account the interlayer overlapping of the n and a bonds yields enhanced values 
of the intrinsic-SOI; about one order of magnitude larger than in single layer graphene 
( ~ 0.1 meV). [42J We would like to emphasize here that such values still somewhat 
weak, compared with those relatively large strengths, of which reportedly, the Rashba 
parameter A# can acquire (similar to the values attained in III-V semiconductors). 
Nevertheless, for generality, the intrinsic SOI has been incorporated in the present 
derivation of the low energy effective Hamiltonian. This will be helpful when considering 
BLG in the extreme limit, i.e. when the intrinsic-SOI rj is much stronger than the 
Rashba-SOI Xr parameter (77 3> \r). However, we shall concentrate here our discussion 
on the bandstructure and Landau levels of BLG for the case A# 3> 77, the opposite limit 
will be treated elsewhere. 

The last term in Eq.(5) arises if an external magnetic field B is present, affecting 
the energetic of the quasiparticles in the form of a Zeeman interaction, Hz, for a field 
perpendicular to the BLG plane it reads 

H z = A(I®a z ), (8) 

where I is a 4 x 4 unit matrix, 2A = g^sB is the Zeeman splitting energy, g is the 
electron Lande factor, and fiB is the Bohr magneton. Note that even at relatively 
high magnetic fields (B = 10 T), the Zeeman splitting it is still somewhat small 
(A ~ 1.1 me V), whiles it is practically negligible at low fields. Note that the condition, 

r] < A < x R <m: 71, (9) 

typically holds at finite fields (B > 0.1 T). This condition will allow us to work safely 
within the low energy theory and derive an effective Hamiltonian for BLG, including 
the extrinsic (Rashba)-SOI, intrinsic-SOI, and Zeeman effect on the same footing. We 
finally should remark that the total Hamiltonian (jSJ is valid near K symmetry point 
only. For the K' point of the Brillouin zone, the Hamiltonian Hk' = Ey-Hi^E" 1 , with 
Y, y = a y ® J, should be used instead. 

2.1. Low energy bilayer Hamiltonian 

Using Lowdin partitioning theory JHJI35] the full 8x8 Hamiltonian H^ can be projected 
through a canonical transformation [46] into a 4 x 4 low energy effective Hamiltonian H 
in an appropriate basis (see Appendix A). The projected low energy Hamiltonian can 
be further expressed in terms of Kronecker products of 2 x 2 matrices and conveniently 
separated into the sum of the Hamiltonians (keeping terms up to I/71), 

-H = H {o) +H {1) +H (2) + 0(l/j 3 ), (10) 

in which the term independent of the interlayer hopping parameter ji reads, 

H {o) = -a z ® (Ua + Aa z ) -n{a ® a z ), (11) 
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where a z is the z— component of the Pauli matrices and a is the 2x2 unit matrix. The 
dominant contribution to % is described by the Hamiltonian T-L^\ given by 

H a, = _l!fo M 2 W + H!Wo -') 83+ , (12) 

71 \tt o y 7i yvr o y 

where we have defined the operator s± = |(c ± °"z)- Without Rashba-SOI (A# = 0), 
Eq. 012jl decouples to the usual effective BLG Hamiltonian obtained within low energy 
theory in the absence of trigonal warping effects. [13] Such term gives rise to well known 
parabolic spectrum of the massive Dirac quasiparticles in BLG. The second term in T-L^ 1 ' 
is linear in momentum and can be viewed as a renormalization of the Rashba coefficient 
due to the presence of the higher bands. Notice that it scales as the inverse of the 
interlayer hopping energy 7^ 

The remaining terms proportional to I/7 2 in Eg. ([TO]) are compacted into the sum 

^ (2) = ElA (2) , with 



h? 



h 2 



(2) 



(2) 




/? 3 
hf 

The low energy effective Hamiltonian % described in Eq. (TTOj) is valid within the 
energy range e < 71. Notice that it will be fairly sufficient to keep only the leading 
order contribution h\ in H^ given the typical smallness of the ratios X'r/jx, AflA/7 2 , 
and XrI]/^ appearing in T-L^ together with the assumption U < 71. Hence the 
description for the low energy (and momentum) effective Hamiltonian will be given 
by u = U^+H^ + hf ) . 

If we further neglect the Zeeman and the intrinsic SOI (A = rj = 0) 
the effective bilayer Hamiltonian with Rashba-SOI written in the atomic basis 
{i>A n ,i>A 2i ,i>B lt ,i>B u } reduces to 





/ -U + £rfir -f3(^f 


\ 






H = 


—U + ^7T t TT —ian^ 
— (3n 2 ia-rr U — ^itit 1 


-/3(vrt) 2 
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\ -/3vr 2 


U-£irir 1 J 




where we have defined the parameters £ = 2L/7 2 , a = 


27 \r and ft = 


= 7i7 2 


7 = 7/7i- 











(13) 



with 
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2.2. Bilayer graphene spectrum with Rashba effect at zero field 



Without magnetic field (B = 0), ^ = h~k_ = H(k x — ik y ) and n = hk+ 
and the eigenvalues of Eq. ( 113ft are readily determined by (Appendix B) 



n[K>x ~\~ %t\j\. 



e^ s {k) 



4 (U - £k 2 ) 2 + k 2 (Ja 2 + Ak 2 P 2 



set 



(14) 



with k = ^Jk 2 + k 2 Here \x = ± describe the electron/hole branch, whiles s — ± 
characterizes its spin chirality. Therefore, the low-energy spectrum consist of four spin- 
splitted bands, two conduction and two valence bands. For the unbias voltage case 
(U = 0) the spectrum reduces simply to 

e^k) = |fc (Ja 2 + A(5 2 k 2 -sa\. (15) 
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Figure 1. (Color online) Low quasiparticle energy spectrum for bilayer graphene 
with Rashba-SOI. The spin degeneracy of the bands at \r = (a) is lifted for Xr ^ 
(b)-(d). As the strength of the Rashba parameter is increased, the symmetry of the 
bands is broken producing a cone chape for the innermost bands (s=-) at high values 
of A fi - 

We note that in contrast with single layer graphene, (47] in BLG the Rashba-SOI 
it is induced a linear spin-splitting in momentum of the conduction and valence bands, 
in close analogy with the Rashba interaction arisen in two-dimensional electron gases 
in semiconductor heterostructures. In addition we observe that the cyclotron effective 
mass at the Fermi energy m c s = k^de^/dk) turns out to be spin- dependent as long 
as Xr 7^ following the relation, 



m 



W 



m 



4/i^/7mA / 3W 1 + 47mA 



1 + 87mA| - sJl + Aim\ 2 



(16) 
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9 



with \p = (3/ a = 7/2A#, and we have expressed the Fermi wave number in terms of the 
2D carrier density via kp = ^prm. Notice that in the limit case A/? — > oo, the cyclotron 
mass m c s — > /J>rn* = fX'ji/2v F , as one expects for unbiased BLG in the absence of SOI. 
However as the carrier density n — > the cyclotron effective mass does not diverges as 
predicted by tight-binding models. 

The normalized eigenvectors \w^J) of T-L corresponding to the electron and hole 
bands /i = ±, respectively, for the case of spin up (s = +) are written in the four-vector 
form (Appendix B) as, 

/ -ie~ 2i<l 'sm{6/2)x+ ) \ 

e-^cos(9/2) X { + 

-i cos(#/2) 



Itfi) 



1 



i + (xT) 



(m)\2 



V 



ic 



-icj> 



(17) 



sin(0/2) ) 



whiles the normalized eigenvectors for the spin down (s - 
specified by 

ie- 2i 't , cos{6/2)x M \ 
I e~ i *mn{6/2)x- ) 



-) electron/hole bands are 



(*0\ 



ICO 



(m)\2 



i + on 



i sin(0/2) 
\ ie~^cos(0/2) 

in which we have defined the dimensionless parameter 



7 



y (m) 

A s 



u-ik 2 + fiJni + {u- m 



K 



(19) 



7'*' 



with 9 = tan 1 (2(3k/a), and <fi is the azimuthal angle of the in-plane wave vector, 
k = k(cos4>,sm(p). The denominator of Eq. (Tl9|) is explictly, 1Z^ S = /ie° s (/c) = |e^ s (fc)|, 



which implies 1Z +S = TZ- S , and therefore the relation xi- +) X. 



(+)-v(-) 



■1 it is always satisfied. 



Without external bias voltage (U = 0), the parameter x° reduces to ±1 for all k 



(±) 



2.2.1. Spin and valley polarization The expectation value of the valley (charge) 
polarization and spin orientation are defined as (t) ms = (ip^ \t\iP^ s ), and (S)^ s = 
{ipfa \S\tf)fc, ), respectively. Here the valley and spin operators are r = cr ® o~ and 
S = f (o"o <S) c), with cr = (da;, a y , a z ) the vector of Pauli matrices. Using the results of 
Eqs. ( fTTl) and ( fl8l) the components for the charge polarization leads to the expressions 



\'x) [IS 
Vy/lis 



2X 



(/■') 



sin6'cos(2(/)), 

(x^y + i 

sin#sin(20), 



2xi 



(x 

(x 






+ 1 
- 1 

+ 1 : 



(20) 
(21) 
(22) 
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whereas the components of the spin polarization (in units of h/2) satisfy 

(S x ),,s = -ssin6 l sin(0), (23) 

(S y )„ s = +ssm6cos((j)), (24) 

(S z ), s = s— f^-cosO. (25) 

1 + (Xs ) 

As it occurs with the standard Rashba SOI in semiconductors, in BLG the 
orientation of the spin-polarization in the plane is always perpendicular to the direction 
of the momentum, (S) ■ k = 0. We notice also that, in contrast with single layer 
graphene, in BLG the dot product (S) ■ t ^ in general. Interestingly, as long there 
is a bias voltage present (U ^ 0), both the spin and valley polarization have a nonzero 
component out of the BLG plane (along the z— axis). However, the absence of bias 
voltage (U = 0) the amplitude of the charge and spin polarization develops k— dependent 
oscillations. Explicitly, |(r)| = \(S)\ = sm.9 for all \k\ ^ 0, and vanishes at k = 0; in 
close analogy with the known result in single layer graphene. |48] From Eqs.(22)-(24) the 
spin-polarization in the unbias configuration can compactly written as 

(S), s = -. 2S(3 [z x fc), (26) 

that is, (S)fj, s is forced to lay on the BLG plane. Clearly, as the Rashba SOI coefficient 
A/j — > 0, i.e. a — > 0, the magnitude of the spin-polarization reaches it maximum value, 
|(S') AiS | — > 1 as the electron/hole spin is conserved. Finally, we notice that Eq. ( j26l) 
for the spin orientation in unbiased BLG is formally identical to that obtained in single 
layer graphene with Rashba SOI. [37] 

3. Landau levels in bilayer graphene with Rashba SOI 

For a magnetic field B^0 perpendicular to the BLG plane, the operators ir and 7T do 
not commute any more since its components fails to do so, and of course, care has to be 
exercised in their ordering. By making the substitution in Eq. (TT3l) of the momentum 
operators in terms of the Bose operators, ir^ = \/2ha) /Ib and n = \/2ha/ls with 
[a, at\ =1, the effective Hamiltonian in the limit of low bias (£/ <C 71) takes the form 

( U u(a) 2 \ 

U iVa id (a,y 

w(at) 2 -iVa) -U 

\ w(at) 2 -U 



H,=- 



n'2 



(27) 



with the notation V = \/2ha/lB and to = 2h 2 (3/l B . The eigenfunctions of Hl can now 
be written in the form, \ijj n ) = (q~ \n — 2), c£ ~ \n— 1), C3 \n), c\ \n + 1)) T , where 
\n) = ^ n are the usual harmonic oscillator eigenfunctions satisfying a)\n) = \Jn + l|n+l) 
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and a\n) = y/n\n — 1). Consequently one can write the expectation value (ip n \H L \ip n ) = 

(<t>n\Hn\<f>n), with 




-H n = ~7 T n , (28) 



and \4> n ) = (cf \c~2 > c 3 i C T ) T satisfying the normalization condition (4> n \4> n ) = 
1, whiles 

( a jft=vj <> V (2 9 ) 

\ iTy/n u)y/n(ri+l) I 

The eigenvalues of ( 128ft leads to the Landau spectrum given by Eq. (2), which, in 
the absence of a bias gate voltage across the layers ({7 = 0), reads (Appendix C) 



e° n . = JL \JnT 2 + 2n 2 u 2 ± nV^u 4 + 4nu 2 T 2 + T 4 , (30) 

v2 

for n > 2. In order to gain further physical insight of the behavior of the Landau levels 
in BLG, it is illustrative to analyze the limit cases at zero, weak(large) Rashba-SOI 
relative to the magnetic field strength, with and without bias voltage. 

3.1. Approximate solutions for U = 0. 

(i) Zero Rashba-SOI (T = 0). In the vanishing Rashba-SOI strength limit, Eq. (1301 
reduces to s° ± = [i\Jn{n ± 1) u B, with u = eh/m*c, and coincides with the LL 
spectrum (double degenerate in spin) reported in the literature for BLG in the absence 
of Rashba-SOI. The linear response with B stems from the parabolic dispersion laws of 
BLG for this case. 

(ii) Weak Rashba-SOI (T/oj <^. 1). At very weak Rashba-SOI strength values (large 
fields), the LL's still evolve approximately linear with B, but shifted by a small energy 
proportional to A|,, described by 

<m± - H (uoB ± i-J Vn(n ± 1), /i = ± (31) 

with r o = ^J2uo /m* (Xr/vf) and n > 2. Since Tl/Au = X 2 R /q\, only for rather large 
Rashba SOI coefficient (Xr ~ 71) strengths gives rise to significant broken degeneracies 
of the electron/hole LLs as described next in the opposite regime. 

(Hi) Strong Rashba-SOI (T/u 3> 1) Alternatively, in the very strong Rashba-SOI limit 
(small fields), the LL level spectrum is well described by s° n ^ + — fiy/nT = fiT ^nB, 
and s° nil _ ~ ii((j0 2 /T o )^n(n 2 — 1)£> 3 / 2 , with n > 2. The change in the field dependence 
from B 1 / 2 of the v = + chiral states to B 3 ^ 2 for the v — — states is unique in BLG and 
give rise to a multiplicity of LL crossings as shown later. 
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In addition, the energy spectrum with LL level index n = 1 and n = are special 
cases, giving rise to three levels, one at zero energy (e - = 0) and two nondegenerate 
levels for n = 1. In the high field limit, T/X R 3> 1, we get £a M + ~ ^^/2{T 2 /Auj + ui B), 
whiles -&l m - ~ fiT y/B in the weak field regime. 

3.2. Approximate solutions for U ^ 0. 

(?) Zero Rashba-SOI in the u> ^ U limit. In this case the quantum states should follows 

e nitx± ~nU + ^jn(n ± 1) u 2 B 2 , (32) 

clearly, in addition to the aperture of an energy gap of 2U between the spin-degenerate 
electrons/hole LLs, the presence of the gate voltage induces a deviation from the linear 
dependence in B occurring at U = 0, to a parabolic behavior with B instead. This is 
also a known result in the literature. [HI [12]. 

(ii) Weak Rashba-SOI in the regime U ^> u> 3> I\ Here the LLs still behave quadratic 
in B but with a spin- dependent shift linear in B given as 

TT un(n±l) ( 2 „ 2 T 2 B\ , . 

e n ^ ± ~vU+ * y m ] lu 2 B 2 ± -|- J , (33) 

the term proportional to V 2 is responsible for the anticrossings of the fan spectrum of 
the v — ± states. 

(iii) Strong Rashba-SOI in the regime U ^> Y ^> cu. In contrast with the case of C/ = 0, 
at very large Rashba-SOI strengths the LL spectrum follows (to leading order) a linear 
behavior with B for the positive chirality states [y = +), whereas for the negative 
[y = — ) quantum states develops a cubic dependence instead. Explicitly they are given 
by, 

e n ^ + ~nU + ^fnB, e n ,,_ ~ fiU + ^n(n 2 - l)B\ (34) 

such drastic change on the field dependence of the LLs with different spin chirality 
states [y = ±) will enhance dramatically the degree of their spin-splitting and on the 
multiplicity of the level crossing as studied in the next section. 

All of this discussed above holds for n > 2. The cases n = 0, 1 are treated separately 
as in the condition for U — 0. In the limit U ^> u ^> T (small Rashba SOI or large 
fields) e lll+ ~ n(U + u 2 B 2 /U), and e lfl+ ~ fi(U + T 2 B/U) for the large Rashba-SOI (or 
small Rashba-SOI) case ([/>r>w), whiles £i_ = —U for both regimes. The case for 
the zero Landau level index is £q- = U, which also holds in both regimes. 

3.3. Spin-polarization of the Landau levels 

The eigenfunctions of Eq. (|28|) for the n— th Landau level of a given electron(hole) band 
\x in biased BLG are given in Appendix C. From Eqs. A0.18J) and the orthogonality 
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Figure 2. (Color online) Low energy spectrum for biased bilayer graphene with 
Rashba-SOI. Here U — 0.050 eV. The bias voltage induced gap decreases as the Rashba 
parameter Ah increase. 



fUJ 



of the oscillator wave functions £„, it follows that, the valley (T {n) )^ u and (S r(n) ) 
spin polarization lying in the plane of BLG vanishes identically for all Landau levels. 
General expressions for the valley and spin-polarization are provided in Eqs. ( 1C.19I) and 
(IC.20J) . We notice that the valley polarization in the perpendicular direction turns to 
be &;— independent, and that for the unbiased case, it vanishes for all LL. Furthermore, 
the z-th component of the spin polarization of the n— th LL with U = 0, reduces to 



W> 



fU> 



— (cos$. 
2 v 



7 ,_ + cos$„ + ) 
2vu 2 



V4w 2 + 4nw 2 r 2 + r 4 ' 



(35) 



where, as before v — ± denotes the plus/minus n— LL of a given \x — ± electron/hole 
branch. In the limit of high field (S^)^ — > =Fl, a full polarization is reached, and the 
states v = ± coincides with the spin-magnetization signs (=f) of the LLs. If the limit 
B — )• is taken, then the spin-polarization (S { z n) )^± — > ±2{oj1/Y 2 )B. 



4. Band structure properties in BLG with Rashba SOI 

The low quasiparticle energy band structure for unbiased bilayer graphene with Rashba- 
SOI, as predicted by Eq.( TT5|) at zero field, is illustrated in Fig.l for different values of 
\r strength. In the absence of Rashba-SOI (Xr = 0) the well recognized parabolic 
spin-degenerated conduction and valence bands touching at its extrema at k = are 
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Figure 3. (Color online) The Rashba-SOI effectively modulates the gap size of biased 
BLG. Inset (a) shows the spectrum of the bands for the critical Rashba-SOI parameter 
that ensures a closing gap, Xr — 71 /v2 with U = 0.025 eV, 71 = 0.22 eV and 
Xji = 12A . 



depicted in Fig. 1(a). For non-zero A# the spin-degeneracy of the bands is broken 



inducing a /c-linear spin-splitting of width A s (k) = \E± iT — E ± 



ah. For relatively 



weak Rashba-SOI (4/3 2 fc 2 3> a 2 ) the band dispersion follows a parabolic behavior, 
£°(k) ~ fi(P 2 k 2 — ^sak), as shown in Fig. 1(b). On the other hand, if the condition 
4(3 2 k 2 <C a 2 holds, then the relation ( fl5l) evolves to a linear spectrum for the inner 
bands (£Z s (k) ~ fxak for s = — ), and to a A;-cubic spectrum for the outermost bands 
(e° (k) ~ fi(3 2 k 3 /a for s = +) as plotted in Fig. 1(d). Similar drastic changes in the 
bandstructure due to large Rashba-SOI strengths were reported earlier numerically by 
van Gelderen et al. [39] within a tight-binding framework; see for instance the low energy 
bands near the Dirac points of Fig. 1(d) of that reference. 

At the intermediate regime (Fig. 1(c)), the innermost bands interpolates from a 
k- linear behavior for electron/hole momentum very close to the Dirac point, to a k- 
cubic dependence for high momentum. In contrast, the s = + bands seems to be 
well described by the cubic spectrum for all values of k. Such remarkable asymmetry 
behavior of the spectrum of BLG with Rashba-SOI is certainly unique, since it is not seen 
in monolayer graphene neither in semiconductors with the Rashba-type of SOI. These 
peculiar characteristics of the spectrum of BLG would have interesting consequences on 
the electronic and spin-transport properties. 

In Fig. 2 we show the bilayer graphene low energy spectrum for finite bias voltage 
(U = 0.050 eV) at various Rashba coupling strengths (Xr = 0,A o ,4A o and 12A D ). As 
in the unbias case, without Rashba-SOI, the bands are spin degenerated (Fig.2(a)). A 
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Figure 4. Spectrum of the LLs for bilayer graphene with Rashba-SOI within the low 
effective low energy theory. A rather unusual characteristics of the LLs is predicted to 
occur at large values of the Rashba-SOI strength. 



gap of 2\U\ at k = is opened between the conduction and valence bands turning BLG 
into a semiconductor. Moreover a band-bending appears at small wave-numbers (low 
momentum) due interplay with the bias gate voltage U. This is the so-called " Mexican- 
hat-like" shape of the lowest energy bands well reported in the literature. From Eq. (IB.8|) . 
in this regime (Xr = 0) and ka < 1 the bands are reasonably well described by 
Sfj,s(k) — n(U — £k 2 + (/3 2 /2U)k 4: ). For nonzero X R (Figs.2(b)-(d)) the spin-degeneracy 
is lifted. As the Rashba-SOI parameter increase the lowest /highest conduction/valence 
bands becomes more warped and the gap tend to decrease as the lowest conduction 
band E ++ evolves from a Mexican-hat like shape to an inverted one, and vice versa for 



the highest valence band E_ + , see Fig. 3(d) calculated using Eq.( 1A.14|) . The behavior 



of the gaps A g+ = E ++ — E_ + (for the innermost bands) and A g _ = E + _ — E (for 

the outmost bands) as a function of the ratio Xr/ji is plotted in Fig. 3 for different bias 
voltages U. The gap A g+ in BLG closes as the Rashba parameter increases, reaching 
its minimum (zero gap) at Xr = 71/V2, to then gradually and linearly opens again as 
Xr/^i is increase up to 1. For Xr/ji > 1 the gap A g+ remains constant. Inset (a) of 
Fig.3 shows the bandstructure for biased BLG with U = 0.025 eV illustrating the zero 
gap condition. An analogous behavior can be seen in the numerical plot depicted in 
Fig. 7(b) of Ref. [39J. Notice that our analytical low effective modeling for the bilayer 
graphene bandstructure allow us to capture also these anomalous behavior of the lowest 
bands (near the Dirac points) under finite bias and relatively large Rashba coupling. 
In addition it predicts that the closing of the gap occurs provided that Xr = 71/V2, 
regardless of the magnitude of the bias voltage applied. 
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Figure 5. The presence of a bias voltage split in two the fan diagram opening a gap 
of 2U between the electron and hole LLs. 



5. Landau level spectrum in BLG with Rashba-SOI 



The Landau level energy spectrum as a function of magnetic field (up to 12 T) for BLG 
with Rashba-SOI according to formula Eq.([2]) is plotted (from n = to n = 41) in 
Fig.4 and Fig.5 for the unbiased and biased case, respectively. The zero gate voltage 
(U = 0) and without Rashba-SOI case shows a LL fan diagram which is linear with 
B and degenerate in spin (Fig.4(a)), as expected, because of the parabolic behavior of 
the energy bands and since there is no spin-dependent mechanism here to break the 
spin-symmetry When the Rashba-SOI is present the spin- degeneracy is lifted inducing 
multiple crossings of the LLs at the Fermi energy, similar as it occurs in semiconductor 
2DEGs with Rashba SOI, and stems because for sizable Xr strengths the LLs with high 
index and with spin chirality s = + have lower energies than those with spin chirality 
s = — as the field is increased. This is basically an intermediate regime between those 
LLs discussed in section 3.1(ii) and 3.1(m). Notice that for a relatively weak intensity of 
the Rashba parameter (Xr = A ) the LLs behave roughly linear with the field (Fig. 4(b)), 
no matter its spin chirality. 

However, for large values (Xr = 4A ) a drastic and unusual change in the LLs 
spectrum arises (Fig.4(c)); the LLs with spin v = + evolves as B l l 2 whiles for those 
with v = — develops a B 3 / 2 dependence as described in section 3.1(m). Such difference 
in the field dependence effectively squeezes the LLs with v — — to lower energies as 
Xr is increased promoting multiple crossings between the LLs. This surprising result 
suggests a strong spin polarization of the Landau levels induced by a significant increase 
of the Rashba parameter. 

The LLs for the biased U ^ and without Rashba coupling case shows instead a 
parabolic behavior with B, and the fan diagram is split by a gap of 2U (Fig.5(a)) as 
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also described by Eg. (132 p . The linear behavior observed for U = and A# = A of the 
LLs transforms as well to a B 2 dependence for U ^ 0, as predicted also by Eg. (133]) . 
The presence of the Rashba-SOI manifest as well as crossings of the LLs as seen in 
Fig. 5(b). The behavior illustrated in Fig. 5(c) at relatively large Xr is similar to that 
seen in Fig. 4(c) but with an open gap and with almost a linear behavior with field for 
v — +, and roughly a cubic dependence in B for the v — — LLs, se also Eg.f l34|) . Such 
asymmetric response at large Rashba-SOI in BLG contrast radically with the linear 
behavior with field that occurs in single layer graphene[37] at the same regime. 

6. Summary 

We have studied the problem of the influence of the Rashba spin-orbit interaction on 
the bandstructure of biased and unbiased bilayer graphene. Using low energy effective 
theory we have derived a low energy Hamiltonian for bilayer graphene in the presence 
of an external magnetic field and spin-orbit interactions. Analytical formulae for the 
energy spectrum of a graphene bilayer with Rashba spin-orbit interaction are obtained. 
We show that for relatively weak Rashba coupling the spin-degeneracy of the electron 
and hole bands is broken inducing a A;-linear spin-splitting, very similar to that found 
in semiconductors heterostructures. At the intermediate strengths of the Rashba effect, 
the innermost bands interpolates from a fc-linear behavior at small momentum, to a 
/c-cubic dependence at high momentum. In contrast, the outermost bands seems to 
be well described by the cubic spectrum for all values of k. For large values of the 
Rashba coupling there is a remarkable warping behavior of the spectrum near the 
Dirac point. Such behavior is unigue in biased bilayer graphene. It is found that 
the bias-induced gap in bilayer graphene decreases as the Rashba is increased, showing 
a behavior resembling a topological insulator transition phenomena. These peculiar 
characteristics of the spectrum of bilayer graphene with Rashba spin-orbit interaction 
may have important conseguences on its electronic and spin-transport properties. 

We also obtained an analytical expression for the Landau levels and spin- 
polarization in biased bilayer graphene with Rashba effect valid in the low bias regime. 
It is further predicted the appearance of an unexpected asymmetric spin-splitting and 
crossings of the Landau levels due the combined effect between the Rashba interaction 
and the bias voltage. These results suggests significant conseguences on the Shubnikov- 
de Hass oscillations and magnetotransport in bilayer graphene as guantum and spin 
Hall effects, under the presence of sizable Rashba spin-orbit interaction in the range of 
few me Vs. 
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In this appendix we derive the low-energy Hamiltonian for bilayer graphene in the 
presence of both, intrinsic and Rashba spin-orbit interaction, as well as Zeeman effect. 
The low energy Hamiltonian is obtained via Lowdin partitioning, also known as van 
Vleck's perturbation theory in the context of atomic physics. |1U E5j H6] First, it 
is convenient to express the Hamiltonian H^ of Eq.(5) in the new spin-dependent 
atomic basis \W K ) = {tpA n ,ipA 2V ^>B^,ipB ll ,^A^,ipA u ,'^B 2t ,ipB 2 J leading to the 8 x 8 
Hamiltonian 

( ~ v 
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with 5+ = , and 5_ = 1. On the other hand 
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(A.3) 



Notice that H s is nothing but the single-layer graphene Hamiltonian with Rashba- 
SOI in the basis {^A mt ,ipA mv tpB mt ,^B im }- 

Bilayer graphene with RSOI, ISOI and Zeeman effect has (in general) eight non- 
degenerated levels. The levels are given by the eigenvalues of the Hamiltonian "H , 



(A.l) 
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(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 

(A.10) 
(All) 
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The eigenvectors |\& M ) of 1-L (with /j, 
the matrix 
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(A.12) 



with C± = cos(i?±/2), S± = sin(i? ± /2) and tan i? ± = ji/(U ± 77). Because of 
the strength of the parameter 71, the energy levels of the subspace of high energy, 
£° a G {E°,E%,E%,E$}, and the energy levels of the subspace with low energy, E° b G 



co co 

\°ia °ja\ 



fo_co\ <ir - 



jb\ 



{£'3, E%, Eg, Eq} are well separated from each other (i.e. 
\£° a — £° b \ ~ j 71 1 ). The low energy Hamiltonian for bilayer graphene can thus be 
obtained through the unitary transformation % = e tS HkC~ lS , in which the S matrix 
elements are given by 
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^ {El - Eo) (E%, - Eoy 

with Sfj, u = (Svfj)^ and W^ u =< ^/^\W\^/ U > . The low energy matrix elements of the 
effective Hamiltonian up to second order in 1/7! are determined by 



7V = K V + w ** 
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e; - El E°, - El 



0(2), 



with r H. ilv = (Hvfj,)^ and /i,// G {3,4,5,6} and v,v' G= {1,2,7,8}. The effective 
Hamiltonian matrix elements reads, 
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These matrix elements form the desired 4x4 low energy Hamiltonian for bilayer 
graphene 
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(A.13) 



which written in terms of suitable Kronecker products leads to Eq. (TTUj) . Note that Eq. 
(IA.13I) was projected on the basis set {ipA 2 i^A 2t ^' l pB l1 -,' l pB u }- 

Energy dispersion. Consider vanishing intrinsic SOI and no magnetic field (n = 
A = 0). To a good approximation we can neglect the off-diagonal terms that goes as 
I/7 2 in the effective Hamiltonian. In such a case the eigenvalues of (1A.13J) are readily 
obtained 
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(A.14) 



here A indicates the electron (+) and hole (— ) branches, whereas s = ± labels 
the spin state chirality, A = A - 4£U, B = (3 2 + £ 2 , A = a 2 + 2fp and T = 
4a 2 /3 2 k 6 + (A 2 + 4p 2 /3 2 ) k A - 2pA(2U - p)k 2 + p 2 (2U - p) 2 , where we have introduced 
the parameters 
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In the limit p ->■ Eq. (1A.14I) reduces to Eq.(14) with A = ji. 



(A.15) 



Appendix B. Eigenvalues of the low energy Hamiltonian 

By squaring the low energy Hamiltonian ( fT3l) . a straightforward diagonalizable system 
is obtained at zero magnetic field (B — 0). 
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with k± = k x ± k y and the eigensystem H? k \x) — £ \\x) yields the eigenvalues 

el ± = {U - £k 2 f + -k 2 (a 2 + 2(3 2 k 2 ± a^a 2 + A(3 2 kA , (B.2) 

whereas its eigenvectors \xj) arranged as column vectors form a unitary matrix 

( -i sin {9/2) i cos (0/2) \ 

e^cos(0/2) e^sin(0/2) 

-i cos {9/2) i sin (0/2) 

y e^sin(0/2) e^cos (0/2) J 

with tan0 = 2/3/c/a, e 1 ^ = {k x + ik y )/k and U^ = U _1 . In the basis of the eigenvectors 
of "Hf,, the 4x4 Hamiltonian Hk is conveniently transformed as follows 
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The eigenvalues of Hk are determined from e^ s {k) = /iW {U — ^k 2 ) 2 + \r s \ , with 
/x = ± for the electron/hole band and s = ± labeling the spin chirality state. Explicitly 
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The eigenvectors of Hk are given by 
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The eigenvectors of Hk are finally given by \ipk) = U|V>jt) leading to Eqs. (TT7|) and 
(USD with x { ^ as given by Eq.flE]). 
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Here we outline the derivation of the Landau levels in biased bilayer graphene 
with Rashba SOI. We follow the same approach used in Appendix B. Squaring the 
Hamiltonian ( |28|) gives the block-diagonal matrix 
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the eigensystem H 2 |<£> n ) = £ 2 n \^Pn) leads to the eigenvalues 

£ l,± = 2 ( 2 C/2 + 2 ™ 2u;2 + rar2 ± n v/4a;4 + 4raw2r2 + r 4 ) , (C2) 

and its corresponding eigenvectors |y?nj) written as column vectors form the matrix, 



V 



( i cos d n+ 

sin -d n+ 

icosfir,- — 2sin$ r 



— % sin -# n+ \ 



V 



sin$ T 







COS$ r 



COS 1?^ 







tan(tf„±) = 2 ^ 2±r2 (C.3) 



/ 



and V'V =1. In the basis of the eigenvectors of % 2 , the 4x4 Hamiltonian 7i n in 
Eq. (128]) is now transformed as follows 

/ U Q n - u n \ 

Q n - -U v n 
v n U Q n+ 
\ u n Q n+ -U J 



u n = v^v 



and 



,)=V- 1 |0n) ; 



(C.4) 



where 

Q n ± = —us [ \Jn(n ± 1) cos f> + cos $_ + \Jn{n =)= 1) sin $ + sin #_ j =p y 7 ^ T cos f9± sin ^ 

(C.5) 



-4u;T- 



n n' 



-1W. 



where we have used the useful relationships 



n > 1 



cos $ nH 
sin $ n+ 
cos t? n _ 
sin $ n _ 



2w 2 


+ t 2 + VK 


2Vn- 


LwT 




Vv 1+ 


v 7 ^ 


— ) 


2^/1 


-i«r 


2/3 2 + r 2 - 


VK 


V 


v 7 ^ 




2uj 2 


-r 2 + v ^; 


2Vn + l 


ooT 




V^>2+ 


v 7 ^ 


— ) 


2y/n 


+ iur 


2/3 2 - r 2 - 


VK 



Vv 



2+ 



v^T 



(C.6) 

(C.7) 
(C.8) 
(C.9) 
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(CIO) 
(C.ll) 

(C.12) 
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with the definitions of the parameters 

M n =4w 2 (u; 2 + nr 2 ) + r 4 , 

2?i± = (An - i)w 2 r 2 + hco 2 + r 2 ± VA4 
v 2± = {An + i)w 2 r 2 + hu 2 - r 2 ± t/K 

Notice that 

u n = —oo ( \Jn(n + 1) cos$ n+ sin#„_ — \Jn(n — 1) sin$ n+ cosf9 n _ J— -v/nrcosi? n+ cos$ n _ = 0, (C.13) 

v n = —u ( \Jn(n + 1) cos$ n _ sint? n+ — \/n(n — 1) sin$„_ cosf9„ + J+A/nrsin , (9 n+ sin$ n _ = 0, (C.14) 

The Landau levels of BLG with Rashba SOI are thus determinated by the 
eigenvalues of H n , which yields 



£ n ,»v = H\JU 2 + \Q nv \ , (C.15) 

for n > 2, being n the Landau level index with v = ± (plus/minus) state of the /x = ± 
electron (hole) band. Using Eqs. ( \C.5\i along with (10.1Q[) - fl0.12)) . the formula ($2§ readily 
follows. The eigenvectors of lri n are given by 

/ 

I0n3> = 



\<f>nl) 



( ° \ 

cos?7_ 

sin?7_ 




V 



\M 



/ COS</ 







V 



V 



\ 

- cos 77+ 

sin 77+ 





JnA 



) 



I — COS ( 

sin (p. 





where the angles r]± and <f>± satisfy, 

' I Qn+ I 



tan 77-1- 



P4 



17 ± ^/U 2 + Q 



and 



tan( 



V 



,(C16) 



M 



± 



U±y/lP + Ql 



,(C17) 



with M_M + = P+P- = — 1. The eigenvectors of H n are thus given by \4> n j) 
V|0 n j), and consequently, the eigenvectors of H n are finally determined by \ip n ) 






in- 



(2)t 

n7 - S»ra-1; 



"./ 






Sn.) 



(4)t 



"J 



The normalized eigenvectors for the states (±) 



of the n— th Landau level of a given electron (hole) band \x are thus explicitly specified 
by 



w:d 



( -isini? n+ sin?7 M £ n _2 \ 

COStfn+sillTfcfn-i 

z sin i? n _ cos 77^ £ n 

y - COS 1? n _ COS 7/>£ n+ i y 



l«> 



/ z cos t9 n+ sin r/ M £„_ 2 \ 

sin^+sin^^i 

-i cos tf n _ cos ip^ £ n 

y - sin t? n _ cos </? M Cn+i y 



(C.18) 



From Eqs. (1C.18I) and the orthogonality of the oscillator wave functions £ n , it follows that 
the components in the plane of BLG of both, the valley and spin polarization vanishes 



» 



(n)\ 



,(n) v 



,(n) v 



identically for all Landau levels, (r^ ) MW = (r^ )^ u = 0, and (Si )^ u = (S x ), 
The valley polarization in the perpendicular direction (along the z— axis) reads 



(rW) 



- cos(2rjf,) , 
sm 7/> - cos z <^ M 



for f = + 
for 1/ = — 



0. 



(C.19) 
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Note that it is &;— independent and that in the limit case of zero bias voltage (U = 0) 
results in rj± = <fi± = f, and hence a zero valley polarization. The z-th component of 
the spin polarization gives on the other hand, 

/Q (u)\ J -cos$„,_cos 2 r^-costf„ + sin 2 ?7 M , u = + , . 

cos i/ n _ cos (pp + cos ij n+ sin r^ M , z/ = — 

which in the absence of bias voltage simplifies to Eq. (|35|) . 
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